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production at HERA. The results presented include some outstanding issues in charm 
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PACS numbers: 12.38.-t, 12.38.qk, 13.60.-r, 14.65.Dw, 14.65.Fy 
1. Introduction 
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^ ' The production of lieavy flavours (cliarm and beauty) in high energy ep collisions at 

Y^' HERA provides new opportunities to study the dynamics of perturbative QCD and to 

^' extract information on the proton and photon structure. 

^ ■ The heavy-flavour cross section at HERA is dominated by photoproduction events 

where the electron is scattered by a small angle, producing photons of almost zero 

^ ■ virtuality. The dynamics of charm and beauty photoproduction can be probed at 

c^ ' HERA in a wide kinematical region, the available centre-of-mass energy being more 

than one order of magnitude higher than in flxed-target experiments. Resolved-photon 
interactions, where the photon fluctuates into parton constituents, which undergo hard 
scattering with the partons from the proton, contribute signiflcantly at HERA energies 
and allow one to study the photon structure. Because of the high cross sections in 
photoproduction, some sort of tagging must be done to allow the data to be recorded 
to tape in the experiment. A charm meson may be tagged and a jet flnder may also 
be used, with a certain number of high transverse energy jets demanded. In the flrst 
case information on the charm meson is available. In the second case the observables of 
the meson, the jet and the correlation between the two may be examined. This sheer 
wealth of variables and currently available statistics is what makes photoproduction most 
interesting. All the previous studies from normal photoproduction may be repeated with 
a charm tag, but thought should be given to original analyses. 

* Summary of the working group on 'Heavy Flavour Physics' at the 3rd UK Phenomenology Workshop 
on HERA Physics, Durham, UK, 20 - 25 September 1998. To appear in the proceedings. 



As a complement, heavy-flavour production at HERA can be studied in deep 
inelastic scattering (DIS). The finding of a scattered positron and the distribution of 
these events in bins of the parton momentum fraction x and the photon virtuality 
Q^ is well documented in the literature. With the additional requirement of a tagged 
charm meson (typically D*) a way is found to measure the charm tagged content of 
the proton structure function F2 {F2). Since -F| is evolved from the gluon content, 
it provides a measurement of the gluon density in the proton. Large statistics in the 
normal F2 measurement and the added measurement of F2 constrain the theory greatly. 
Inclusive F2 is predicted to contain between 10% and 25% charm and it must therefore 
be properly taken into account in global analyses of structure function data. Differential 
cross sections for D* production in DIS have also become available and allow a more 
detailed study of the charm production mechanisms than is possible with the inclusive 
F2 measurement alone. 

The production of heavy-quark bound states has been the subject of intense study 
during the past few years. Exciting phenomenological developments followed from the 
application of non-relativistic QCD (NRQCD), an effective theory that disentangles 
physics on the scale of the heavy-quark mass, relevant to the production of a heavy- 
quark pair, from physics on the scale of the bound state's binding energy, relevant to 
the formation of the quarkonium. The NRQCD approach implies that so-called colour- 
octet processes, in which the heavy-quark-antiquark pair is produced at short distances 
in a colour-octet state and subsequently evolves non-perturbatively into a physical 
quarkonium, must contribute to the cross section. Although gluon fragmentation into 
colour-octet quark pairs appears as the most plausible explanation of the large direct 
ip production cross section observed at the Tevatron, NRQCD factorization is still not 
definitely established on a quantitative level. The analyses of inclusive charmonium 
production at HERA offer unique possibilities to test general features of the quarkonium 
production mechanisms and to establish the phenomenological significance of colour- 
octet processes. 

In the following we will discuss various aspects of charm and beauty production at 
HERA, organized into the three above-mentioned areas: photoproduction, deep inelastic 
scattering and quarkonium production. Additional information may be found in fl]-^. 

2. Photoproduction 

The photon structure is being probed in new kinematic regimes at HERA. Tagging of 
charmed mesons or semileptonic decays in photoproduction events is providing a wealth 
of information on charm in the photon and the processes involved in its production. 
Thanks to the overall higher cross section for photoproduction events, there exists 
a large statistics sample from the 1996 and 1997 data sets alone (37 pb~^). Since 
photoproduction of charm and beauty is much less affected by non-perturbative and 
higher-order effects than hadropro duct ion cross sections, the HERA measurements can 
provide important tests of the heavy-flavour production dynamics. 



The photon-proton cross section may be written as a sum of a direct and a resolved- 
photon contribution: 

(^ipiQ'^) = Yl J^^pfi'i^P^ fip)ai^{xp, as{fJ^R),fip, fJ^R, fJ^-y, Q"^) 

i 

+ J2j dxpdx^ff{xp, fi'p)fj{x^, n^)aij{xp, x^, asifip), fi'p, h'r, h^, Q^), (1) 

where ff{xp,fj,p) and fj{x^,fi^) denote the density functions for parton i and j in the 
proton and in the photon, respectively. The renormalization scale is labelled fiR, and 
Q^ is the hard scale of the interaction. The factorization scales /ip and fi^ are related 
to the subtraction of divergences associated with the collinear emission of partons from 
the incoming parton in the proton and in the incoming photon, respectively. 

The factorization of collinear initial-state singularities has to be performed at NLO 
where Feynman diagrams like Figure ^ contribute. Although this diagram forms part 



Figure 1. A NLO direct photoproduction process. 

of the 0{aal) cross section, if the virtuality k^ of the propagator labelled in the figure 
is less than /i^, then its contribution is already included in the LO-resolved process. 
Hence the separation of direct and resolved reactions is ambiguous beyond LO and 
only the sum is meaningful. This is explicitly shown in (|l]), where the factorization 
scale /i^ appears in both the direct and the resolved contributions. Total cross sections, 
single-inclusive distributions as well as heavy-quark correlations have been calculated 
to next-to- leading order accuracy 0]. 

In higher orders, potentially large terms ~ asln(py/?7i^) arise from the collinear 
branching of gluons or photons into heavy-quark pairs or from the collinear emission 
of gluons by a heavy quark at large transverse momentum. For pt ^ rric these 
terms may spoil the convergence of the perturbation series and should be resummed 
by absorbing the collinear In (p^/?7i^) -singularities into heavy-quark parton densities 
and fragmentation functions. Applications of this approach to charm and D* 
photoproduction at pt > rn^ can be found in |^. For a comparison of the different 
theoretical schemes, see 0. 

Experimentally, the first topic to be studied, in general, is the production of 
charmed mesons, which may be identified by their decay products. Both HI and ZEUS 
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Figure 2. The pseudorapidity distribution of the tagged D* for p^ > 2, 3, 4, 6 
GeV. The points are ZEUS data for the 1996 and 1997 data sample 0. The 
inner error bars are the statistical errors and the outer error bars are the 
statistical and systematic errors added in quadrature. The curves are different 
theoretical NLO predictions yj5|. The variable e is the non-perturbative 
parameter in the Peterson et al. fragmentation function |lO|1 . 



have shown results on the D* meson in photopro duct ion ^^. This channel is favoured 
because of its clear signal. A comparison of ZEUS data and different NLO theory 
predictions for the D* pseudorapidity distribution is shown in Figure 0. The data lie 
above the NLO theory when the D* meson travels in the forward direction, as is also 
seen in the DIS case. Although the discrepancy might not be considered very significant 
at present, some interesting suggestions were made for its possible causes. They are 
listed here and discussed below: 



m 



(i) 'String effects' between the proton remnant and the c quark |[Tl|l . 

ii) Low-x-y non-perturbative gluon in the photon |jl2|. This would give rise to more 
charm and is currently poorly constrained by LEP data. 

Effects attributed to the Peterson et al. fragmentation function used. 

(iv) More beauty than predicted by the theory. 

Item (i) can be examined by comparing the available LO Monte Carlo programs to 
the data. In the fragmentation process a string would be linked between the charm 
quark and the forward travelling proton remnant. The tension in this string could be 
responsible for pulling the quark (and hence the meson) forward. This effect is not 
implemented in the NLO calculations, so LO Monte Carlo would be more forward than 
NLO calculations. This comparison was carried out at the recent HERA Monte Carlo 
Workshop |jl3| and a better agreement is seen between the shape of the LO curves and 
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Figure 3. The pseudorapidity distribution of the tagged D* for p^* > 2 GeV 
and pj^ > 4 GeV. The points are ZEUS preHminary data for the 1996 and 
1997 data sample. The inner error bars are the statistical errors and the outer 
error bars are the statistical and systematic errors added in quadrature. The 



sohd histogram is HERWIG flj] and the dashed histogram is PYTHIA [|l| 



both having been area normalised to the data. 

the data (see Figure ^. Hov^^ever, there are many other effects implemented in the 
Monte Carlo (including item (iii)). The process of separating them is still ongoing. 



Previously ZEUS demonstrated [ll^] its sensitivity to the low-x^ regions, by the use 
of the observable 



_ Z^jets ^T ^ 



MBS ^ -Jets^, ^ ^2) 



where the sum is over the two jets with highest transverse energy. In LO, x°^^ can 
be interpreted as an observable related to the fraction of the photon's momentum 
participating in the hard-scale interaction. At higher x^^^ (the usual cut is 0.75) the 
direct processes dominate, and at lower x^^^ the resolved processes dominate. If a 
charm meson or semileptonic decay is then tagged in these events, statements can be 
made about the charm component of the photon. This is of interest because charm in 
the photon may come from the low-x-y non-perturbative gluon, a parton density of the 
photon still badly constrained by LEP PT| . 

Demanding a jet also provides another, even harder, scale to the calculation. 
Typical experimental jet transverse energies are of the order of 6 or 7 GeV, well above 
the charm mass. When Q^ ^ rn^, charm should behave as a light quark. At the scales 
probed with present data it is unclear to what extent this is true, and the problem arises 
of how to treat the intermediate Q^ region. Similar theoretical issues arise for charm 
production in DIS and will be discussed in the next section. 
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Figure 4. The x^^"^ distribution for tagged D* decays from charm, beauty and 
hght quarks. The points are ZEUS prehminary data points for the 1996 and 
1997 data sample. The inner error bars are the statistical errors and the outer 
error bars are the statistical and systematic errors added in quadrature. In 
the top plot the shaded histograms show the direct and resolved contributions 
as predicted by HERWIG 5.9. The solid line is the two contributions added 
together. In the bottom plot the same data is shown (rebinned into four bins) 
against a NLO calculation 0] with two choices of parameters. 



ZEUS have published a measurement of x. 
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with a tagged D* meson, shown 



in Figure ^. This shows a clear tail to low a;^^^, which in terms of LO Monte Carlo 
programs can only be explained by a large (40%) component of LO-resolved processes. 
This is heavily dependent on the amount of charm in the parton density function used. A 
NLO calculation, performed in a factorization scheme where heavy quarks are exclusively 
generated at short distances from light-parton scattering, cannot describe this long tail 
and falls short of both the data and the LO Monte Carlo. This could be due to a large 
non-perturbative part of the gluon in the photon, or could be attributable to effects 
from the fragmentation function used. 

Heavy-hadron fragmentation functions were discussed as an interesting area of 
study and not just as a possible explanation as item (iii) on the above list. The extraction 
of fragmentation functions from experimental data and the value of the corresponding 
non-perturbative parameters are closely related to the description of the perturbative 
part of the production cross section. Different theoretical schemes for calculating heavy- 
quark production require different heavy-hadron fragmentation functions. (See [|I^ for 



recent NLO extractions of D* fragmentation functions from e+e" data.) Suggestions 
have been made for defining and directly measuring fragmentation functions at HERA 
||19|| . However, the invariant mass of the subprocess cannot be determined directly 
and has to be reconstructed from observables in the final state (e.g. jet variables). It 
therefore makes sense to define the fragmentation parameter in terms of the jets. One 
example would be: 

P-p(jet) (i?-p.)jet' ^' 

where P is the four momentum of the proton and p{D*) and p(jet) are those of the meson 

and the jet, respectively. A modified fragmentation function may then be defined as 

1 da^"^-'^' (z) 
D-,.^o-(z) = ^.^.. ,/ ' ■ (4) 

Only a convolution of the c —^ D* and c -^ jet fragmentation functions may be 
measured, so direct comparison with the LEP data is impossible. With suitably large 
statistics, it may be possible to do a search for scaling violations with HERA data 
alone. Work is now proceeding on the measurement of charm fragmentation functions 
at HERA. 

Finally, the subject of open beauty production was brought up. HI already had 
released their preliminary measurement and comparison with the LO theory [^]. HI 
quote, for the kinematic range Q^ < 1 GeV^, 0.1 < y < 0.8, px > 2.0 GeV and 
35° <e^' < 130°, the preliminary visible cross section of aj'^_^^^^^^^ = 0.93±0.08l[j:[]^^ nb 



and also quote a predicted cross section from AROMA 2.2 |21|] of o'^^^^^,g;,j^ = 0.191 nb. 



This states that the preliminary measurement of the cross section of open beauty was 
up to a factor of 5 higher than the LO theory as predicted by the Monte Carlo AROMA. 
It should be noted that AROMA is designed to produce only a direct component. The 
magnitude of the resolved component is not taken into account in the comparison. 

At the time of the workshop, no result from ZEUS was forthcoming. However, 
before going to press, a preliminary result was released of a beauty measurement by the 



study of semileptonic decays to electrons |^. The measurement is an inclusive electron 
measurement 

e~^p — >• e~ + dijets + X 

in the kinematic region: 

• Q2 < 1 GeV^ 0.2 <y< 0.8. 

• Two jets with Et > 7, 6 GeV and |?7| < 2.4. 

• An electron in the final state with px > 1.6 GeV and \i]\ < 1.1. 

This measurement is inclusive of all beauty, charm and light-quark decays, but with 
all detector effects removed (e.g. conversion electrons). The x^^^ distribution (Figure |^) 
shows, as with the D* study, a clear tail to low x°^^, which needs a significant fraction 
of resolved LO Monte Carlo to describe it. The x^^^ distribution also peaks at high 
values, which is consistent with the observation of direct processes. The preliminary fit 
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Figure 5. The x^^^ distribution for tagged semileptonic decays from charm, 
beauty and hght quarks. The points are ZEUS prehminary 1996 and 1997 
data points. The inner error bars are the statistical errors and the outer error 
bars are the statistical and systematic errors added in quadrature. The shaded 
histograms show the direct and resolved contributions as predicted by HERWIG 
5.9. The solid line is the two contributions added together. 



on the fraction of the LO resolved contribution yields 35 ± 6% and agrees well with the 
HERWIG 5.9 [Q prediction of 40%. The agreement in shape between the data and the 
LO Monte Carlo is good. 

Next the p^"^ distribution was examined. In this study p^^' is the momentum of the 
electron transverse to the jet axis of the closest jet (the jets were reconstructed using a 
clustering algorithm, KTCLUS [^). The electrons from beauty decays are expected to 
dominate at high p^' because of the larger beauty mass. The obtained distribution is 
shown in Figure |^. A fit was done allowing the fraction of beauty to vary with respect 
to the charm and light quarks contribution (added in the ratios given by HERWIG). 
A preliminary beauty fraction of 20 ± 6^]'^% is needed to fit the data (using the 40% 
resolved component), in agreement with the HERWIG prediction of 17%. 

ZEUS quotes a preliminary visible cross section for beauty production of 



^^(^ 



P 



e~ + dijet + X) 



39 ±111?^ pb 



in the kinematic region described above. 

Finally, a preliminary overall normalization factor of 3.7 is necessary to fit the 
Monte Carlo to the data. This normalization factor is for the complete sample and not 
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Figure 6. The p^? distribution for tagged semileptonic decays from charm, 
beauty and Ught quarks. The points are ZEUS prehminary 1996 and 1997 
data points. The inner error bars are the statistical errors and the outer error 
bars are the statistical and systematic errors added in quadrature. The shaded 
histograms show the beauty and charm added to other components as predicted 
by HERWIG 5.9. The solid line is the two contributions added together. 



just for the beauty component. It should be noted that these results depend on the 
parton density functions used in the analysis (for the photon GRV-LO p4| and for the 
proton CTEQ4 [^) and also on the charm and beauty masses used {rric = 1.55 GeV 
and rrib = 4.95 GeV). It should also be noted that for normal dijet photoproduction a 
factor of 1.8 is generally needed to normalize the Monte Carlo to the data. In the D* 
with dijets study, a factor of 2.6 was necessary. 

There is clear consistency between HERWIG 5.9 and the ZEUS data in the 
percentages and p^" '"distribution shapes. However, they disagree with the large 
normalization factor. Taking all the caveats in the previous paragraph into account 
and the different measured kinematic regimes, it is hard to make a statement about the 
agreement between the experiments at present. A comparison of the visible cross section 
with a NLO calculation as well as alternative experimental techniques using microvertex 
detection should allow a clarification of the issue in the future. 



10 
3. Charm Production in DIS 

The production of charm quarks in DIS has become an important theoretical and 
phenomenological issue. The charm contribution to the total structure function F2 
at small x, at HERA, is sizeable, up to ~ 25%; through the charm contribution to 
the scaling violations the treatment of charm also has a significant impact on the 
interpretation of the NMC data. Besides being an interesting theoretical problem due 
to the presence of two hard scales, ttIc and Q, a proper description of charm-quark 
production in DIS is thus required for a global analysis of structure function data and 
a precise extraction of the parton densities in the proton. 

At scales Q ^ rric, charm production in DIS is calculated in the so-called fixed 
flavour number scheme (FFNS) from hard processes initiated by light quarks (m, d, s) and 
gluons, where all effects of the charm quark are contained in the perturbative coefficient 
functions. The FFNS incorporates the correct threshold behaviour, but for large scales, 
Q ^ mci the coefficient functions in the FFNS at higher orders in as contain potentially 
large logarithms \n^{Q'^ /ml), which may need to be resummed. Such a resummation can 
be achieved by including the heavy quark as an active parton in the proton. The simplest 
approach incorporating this idea is the so-called zero mass variable flavour number 
scheme (ZM-VFNS), where heavy quarks are treated as infinitely massive below some 
scale Q ~ mH and massless above this threshold. This scheme has been used in global 
fits for many years, but it has an error of 0{m\/Q'^) and is not suited for quantitative 
analyses unless Q ^ mn- Considerable effort has been made to devise a scheme for 
heavy-flavour production that interpolates between the FFNS close to threshold and the 
ZM-VFNS at large Q. Here, we will focus on two such schemes, the Aivaziz-Collins- 
Olness-Tung (ACOT) |2^ and the Thorne-Roberts (TR) [^ scheme, which have been 
used in recent global analyses of parton distributions.^ These generalized VFNSs include 
the heavy quark as an active parton flavour and involve matching between the FFNS 
with three active flavours and a four-flavour scheme with non-zero heavy-quark mass. 
They employ the fact that the mass singularities associated with the heavy-quark mass 
can be resummed into the parton distributions without taking the limit tuh — > in 
the short- distance coefficient functions, as done in the ZM-VFNS. To all orders (and 



neglecting intrinsic heavy-quark contributions ||28[) the FFNS, the ACOT scheme, and 
the TR scheme are identical, but the way of ordering the perturbative expansion is not 
unique and the results differ at finite order in perturbation theory. A NLO comparison 
of the predictions for F| within the FFNS [g9|,|g, the ZM-VFNS and the TR scheme 
is shown in Figure |^. While the ZM-VFNS does not provide an adequate theoretical 
description at low values of Q, the difference between the FFNS and the TR-scheme is 
moderate, typically ^ 20%. A similar behaviour is observed for the ACOT-VFNS. 

Experimentally, ZEUS have released their F2 measurement using their full 1996 



and 1997 data samples |^I|. HI have released their measurement using only the 1994 
to 1996 data sample ^. Figure | shows the Q^ and x dependence of F2 compared to 

t For theoretical details and the discussion of other schemes, see B. 
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Figure 7. Charm quark structure function, F^{x,Q'^) for x = 0.05 calculated 
in the TR-VFNS, FFNS and ZM-VFNS. See i% for details. 



the NLO FFNS calculation ||30[. The measurements are in good agreement with the 
theoretical prediction. Given the experimental error and the theoretical uncertainty 
due to the variation of the charm-quark mass, as shown in Figure §, present data on 
F2 cannot discriminate between the FFNS and the ACOT and TR schemes. Global 
structure function fits, however, seem to favour the VFNSs over the FFNS [^. 

The resummation of ln\Q'^/ml) terms at higher orders in the FFNS may not 
be necessary for the present analysis of neutral-current structure functions, but it is 
required for a stable theoretical prediction of the charged-current structure function 
-^3^ [i3- ^s analysed in |^, the higher-order corrections to the lowest-order flavour 
excitation mechanism W^ + s —>■ c/c are large and dominated by VT-boson-gluon fusion 
W^ + g ^ c/c + s/s at X ^ 0.1. The extraction of the strange quark density from the 
charged-current structure functions at small x thus requires a good knowledge of the 
gluon density.|] 

Another source of potentially large higher-order corrections to charm production 
in DIS comes from the emission of soft gluons in the threshold region. This effect 
has been studied within the FFNS for charm production through photon-gluon fusion 
The fact that heavy-quark production at HERA energies may be sensitive 



m 



35 



to threshold effects is due to the large gluon density at small x, which enhances the 
contribution of soft-gluon emission near the elastic limit. An all-order resummation of 



threshold logarithms has been performed in |^5[ at next-to-leading logarithmic accuracy; 
an approximate NNLO result has been derived for Fg and the single-particle inclusive 
distribution dFg /dp^- Figure ^ shows the effect of the approximate NNLO corrections 
for the inclusive structure function. Plotted are the i^-factors F2^(NLO)/F2'^(LO) (solid 

t For a discussion of charged-current charm clcctroproduction in the ACOT scheme, see [p4[. 
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Figure 8. The dependence of i^2 o^ Q^ (left) and x (right) compared to the 



NLO FFNS calculation |3C] using the GRV94 parton densities [gl. 



line) and F|(NNLO)/F|(NLO) (dashed line). It can be concluded that the soft-gluon 
effects are well under control for x <^ 0.01, where the size of the NNLO corrections, 
evaluated at the central scale fi = rric, is negligible. It has also been shown in |3^ that 
the inclusion of NNLO terms stabilizes the theoretical prediction for Fg by reducing the 
scale dependence of the NLO result. 

Charm production in DIS has not only been studied for inclusive structure functions 
but differential cross sections for D* production are available as well P,pT[]. Differential 
distributions allow a study of the underlying production mechanism more detailed than 
what is possible with the total cross section alone. A comparison of recent HERA data 
for differential D* cross sections in DIS with the NLO FFNS calculation PDi shows good 
overall agreement, with a possible slight excess of data in the forward region for large 
f]{D*) and in the central x{D*) distribution (see Figure |10|) {x{D*) is the fractional 
momentum of the D* in the 7*p frame). We have verified that the excess is not driven 
by the small-Q region of the data, so a possible resolved photon contribution can be 
ruled out. Hadronization effects may well be responsible for the deviation from the 
parton-level calculation, as discussed for D* photoproduction in the previous section. 

Recently, new theoretical results for differential distributions in the ACOT scheme 
36[| . An event generator Monte Carlo program has been 



have become available 

constructed, which allows to produce arbitrary differential distributions and to impose 



experimental cuts. Figure 11 shows a comparison between the ACOT prediction for 
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Figure 9. The x-dependence of the iT-factors Fg (NLO)/F|^(LO) (soUd hne) 
and F2''(NNL0)/F|(NL0) (dashed line). NNLO result in the approximation of 
35i. 



the pt distribution in D* DIS and recent HI data 0. Given the uncertainty in the 
choice of the non-perturbative D* fragmentation function, which affects the pt slope 
of the theory, the agreement is good. In Figure |ll] the Peterson et al. form [jlOl has 
been adopted, with e = 0.078. A smaller value for e, as suggested by recent analyses 
of e^e~ data [^, would increase the cross section at large pt and thereby improve the 
theoretical description of the data further. A more detailed comparison between HERA 
data and VFNS predictions will be carried out in the future. 

Experimentally much effort is now being put into reducing the systematic and 
statistical errors, to gain as much as possible from the current data samples. First, 
detailed systematic studies are in progress for standard F2, and these will eventually be 
applied to the measurement of F2. Some of the systematic error on F2 is because of 
the extrapolation, due to the forward excess seen by both collaborations BJ3T[]. Work 
is under way to implement 'string effects' described in the previous section into NLO 
calculations. These appear to correct the calculation so as to agree with the data in the 
forward direction. 

The statistical errors may not be much improved by 1998 and 1999 data taking, 
so work is under way to examine different channels. Already ZEUS have shown a 
tagged semileptonic decay channel |3^, and other channels are under investigation (e.g. 
D* -^ KmnnTs). Since these are statistically independent samples, it may be possible 
to combine them to reduce the statistical errors. Finally it is hoped that the new HERA 
luminosity upgrade planned for next year and the inclusion of results from microvertex 
detectors into the analyses will provide vastly increased statistics as well. 
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Figure 10. The differential D* cross section measured by ZEUS [pl| with 
respect to Q^, x, the 7*p centre-of-mass energy W, the transverse momentum 
and pseudorapidity of the D* and x{D*). The shaded band shows the 
predictions of [pG] including the effect of varying rric between 1.2 and 1.6 GeV. 



4. Quarkonium Production 

The production of quarkoniuni at liigh-energy colliders has been the subject of 
considerable interest during the past few years. New experimental results from pp, ep 
and e^e^ colliders have become available, some of which revealed dramatic shortcomings 
of earlier quarkonium models (see 



39 



for recent reviews). In theory, progress on 
factorization between perturbative and the quarkonium bound state dynamics has been 
made. The 'colour-singlet model' (CSM) has been superseded by a consistent and 
rigorous approach, based on non-relativistic QCD (NRQCD) |^, an effective field 
theory that includes the so-called colour-octet mechanisms. On the other hand, the 
'colour evaporation' model (GEM) ^1|] of the early days of quarkonium physics has 
been revived [^. Despite these developments, the range of applicabihty of these 
approaches to the practical case of charmonium is still subject to debate, as is the 
quantitative verification of factorization. The problematic aspect is that, because the 
charmonium mass is still not very large with respect to the QCD scale, non-factorizable 
corrections may not be suppressed enough, if the quarkonium is not part of an isolated 
jet, and the expansions in NRQCD may not converge very well. In this situation cross 
checks between various processes, and predictions of observables such as quarkonium 
polarization and differential cross sections, are crucial in order to assess the importance 
of different quarkonium production mechanisms, as well as the limitations of a particular 
theoretical approach. 
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Figure 11. Differential D* DIS cross section. Compared are recent HI data 
[1 with the theoretical prediction in the ACOT-VFNS l3|]. 



According to the NRQCD formalism, inclusive quarkonium production can be 
viewed as a two-step process, where the production of a heavy-quark pair QQ[n\ 
in a certain angular momentum and colour state n is described by a perturbatively 
calculable short-distance cross section dan and the subsequent quarkonium formation 
is parametrized by a non-perturbative matrix element (O^). The magnitude of the 
non-perturbative transition probabilities is determined by the intrinsic velocity v of the 
bound state and the production cross section can be expressed as a double expansion 
in as and v. 

In the CEM and the related 'soft colour interaction' (SCI) model ^^, the cross 
section for a specific charmonium state is given as a fraction of the inclusive cc production 
cross section integrated up to the open charm threshold. The production fraction is 
assumed to be universal in the CEM, while in the SCI model it depends somewhat on 
the partonic state and the possible string configurations. In both models the colour and 
spin quantum numbers of the intermediate cc state are irrelevant, and gluon emission 
during hadronization is assumed to be unsuppressed. 

The analysis of charmonium production in high energy ep collisions at HERA 
appears to be a powerful tool to test the NRQCD approach and to discriminate NRQCD 
from the CEM or SCI model. The production of charmonium at HERA is dominated by 
photon-gluon fusion, where the photon interacts directly with a gluon from the proton. 
Resolved processes are expected to contribute significantly to the lower end-point of 
the J/ip energy spectrum and might be probed at HERA for the first time in the near 
future. Higher-twist phenomena such as diffractive charmonium production |Q are not 
included in the leading-twist calculations of NRQCD and the CEM: it can be eliminated 
by either a cut in the J/ip transverse momentum pr ^ 1-2 GeV, or by a cut in the J/tp 
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energy variable z = pp ■ k^ / Pp ■ k^ <, 0.9 (in the proton rest frame, z is the ratio of the 
J/ip to 7 energy, z = E^/E^). 

Within NRQCD and at leading order in the intrinsic velocity v of the bound state, 
the inelastic photoproduction of J /if) conies from an intermediate cc pair in a colour- 
singlet ^5*1 state and coincides with the colour-singlet model result. (The notation for the 
angular momentum configuration is '^^^^Lj with S, L and J denoting spin, orbital and 
total angular momentum, respectively.) Cross sections |^^, polar [Q, and polar and 



azimuthal |^^ decay angular distributions have been calculated for the direct-photon 



contribution. The angular integrated cross section is known to next-to-leading order in 
as 1^. Relativistic corrections to the colour- singlet channel at 0{v'^) were calculated in 



49| and were found to be small in the inelastic region z <, 0.9, p^ ^ 1 GeV.|§]The colour- 
singlet contribution, including next-to-leading corrections in ag, is known to reproduce 
the current J/ip photoproduction data [^ adequately. But there is still a considerable 



amount of uncertainty in the normalization of the theoretical prediction, which arises 
from the value of the charm quark mass and the wave function at the origin, as well as 
the choice of parton distribution functions and renormalization/factorization scale. 

At order f^ ~ 0.05-0.1 relative to the colour- singlet contribution, the J/ip can 
also be produced through the intermediate colour-octet configurations ^Si, ^Sq and 
^Pj. In the inelastic region, they have been considered in |Q for the direct photon 
contribution and in [Q for resolved photons, in which case the photon participates 
in the hard scattering through its parton content. Colour-octet contributions to the 
total photoproduction cross section (integrated over all z and pr) are known to next-to- 
leading order ^3| . The polarization of inelastically produced J/ip due to these additional 



production mechanisms has been calculated recently 

The colour-octet production channels are kinematically different from the colour- 
singlet one, because the ^5*0 and ^Pj configurations can be produced through t-channel 
exchange of a gluon already at lowest order in as- (For the ^Si octet this is true for 
the resolved process.) This leads to a significantly enhanced amplitude, in particular in 
the large-z region. The colour-octet contributions to J/ip photoproduction are indeed 



strongly peaked at large z |£2[. Such a shape is not supported by the data, which 
at first sight could lead to a rather stringent constraint on the octet matrix elements 
(OsH), n G {^S'c'^Po}, and to an inconsistency with the values obtained for these matrix 
elements from other processes. However, it is premature to interpret the discrepancy 
between the theoretical predictions for colour-octet contributions and the HERA data 
at large z as a failure of the NRQCD approach itself, as long as not all theoretical 
uncertainties that reside in the calculation of the short-distance cross sections have been 



assessed carefully. Potentially large QCD effects include higher- twist contributions |56[1 , 

higher-order corrections and intrinsic transverse momentum of the partons in the proton 

57| . Moreover, the peaked shape of the z-distribution is derived neglecting the energy 



§ This result has been called into question by a recent analysis JS^I, which finds a sizeable 0{v^) 
correction in the large-z region even after a cut in the transverse momentum, pt ^ 1 GeV, has been 
applied. 
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transfer in the non-perturbative transition cc[n\ -^ J/ip + X. In reality the peak may 
be considerably smeared |5^ as a consequence of resumming kinematically enhanced 
higher-order corrections in f ^ and no constraint or inconsistency can be derived from the 
end-point behaviour of the ^-distribution at present. As a consequence, the role of octet 
contributions to the direct process remains unclear. The resolved photon contribution, 
on the other hand, could be entirely colour-octet-dominated |§^§5] . The ^-distribution 
should then begin to rise again at small z, if the colour-octet matrix elements are as 
large as suggested by NRQCD velocity scaling rules ||60|j40[| and fits to hadroproduction 
data. 

The CEM and the SCI model also predict a strong colour-octet enhancement at 
small and large values of z, again due to t-channel gluon-exchange diagrams [^ . As in 
the case of the NRQCD analysis, the impact of QCD effects (higher-twist contributions, 
higher-order corrections and intrinsic transverse momentum) has to be analysed before 
a final conclusion can be drawn. Since the qualitative behaviour of the J/i(j energy 
distribution is similar in all existing approaches that include colour-octet contributions, 
it is not well suited to discriminate between NRQCD, CEM and SCI. 

The most powerful way to test NRQCD against the CEM/SCI is the analysis of 
charmonium polarization. The polarization analysis in NRQCD ||6^ is based on the 
symmetries of the NRQCD Lagrangian, of which spin and rotational symmetry are 
crucial, and the charmonium polarization can be calculated from the spin orientation of 
the intermediate cc state. In contrast to the NRQCD approach, the CEM/SCI predict 
charmonium to be produced always unpolarized. The models assume unsuppressed 
gluon emission from the cc pair during hadronization, which randomizes spin and 
colour. This assumption is wrong in the heavy-quark limit where spin symmetry is 
at work and soft gluon emission does not flip the heavy-quark spin. Nonetheless, 
since the charm-quark mass is not very large with respect to the QCD scale, the 
applicability of heavy-quark spin symmetry to charmonium physics has to be tested by 
confronting the NRQCD predictions to experimental data. Two instructive examples 



for J/^p photoproduction at HERA are shown in Figure |T2| |^, namely the polar- 
angle parameter A and the azimuthal-angle parameter u in the leptonic J/ip decay, 
as a function of the energy fraction z and the transverse momentum px, respectively. 
The polarization signatures predicted from colour-octet processes within NRQCD are 
distinctive and can also be used to extract information on the relative weight of colour- 
singlet and colour-octet contributions. Polarization signatures have also been studied 
for T production in pN collisions at HERA-B [Q. NRQCD predicts the T mesons to 
be produced transversely polarized, again in contrast to the CEM/SCI, which imply 
unpolarized quarkonium. 

Another powerful way of discriminating NRQCD against the CEM is the 
measurement of Xc photoproduction. In NRQCD the a{xc)/<^{J/ip) ratio is process- 
dependent and strongly suppressed in photoproduction with respect to hadroproduction. 
On the other hand, the assumption of a single universal long-distance factor in the 
CEM implies a universal a{xc)/o'{J /^j) ratio, which is not, however, supported by 
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Figure 12. Colour-singlet (CS) and colour-octet (CO) contributions to the 
polar and azimuthal angular parameters in the leptonic J/ij: decay. For details 
and the definition of the Gottfried-Jackson and Collins-Soper reference frames 



see 55 



the comparison of charinonium production in fixed-target hadron collisions and early 
pfiotopro duct ion experiments. A search for Xc production at HERA, resulting in a cross 
section measurement or upper cross section limit would be crucial to settle this issue. 

Recently, first experimental results on inclusive J /ip electroproduction have been 
presented |B^] and compared to a leading-order NRQCD calculation including colour- 
singlet as well as colour-octet contributions |6^. At large photon virtuahties. 



the analysis of J/ifj leptoproduction is under better theoretical control than 
photoproduction, since higher-order corrections and higher-twist contributions should 
become less important as Q^ increases. Also, at sufficiently large Q^, diffractive 
contributions are expected to be suppressed. The experimental data for differential cross 
sections in the kinematic region Q^ > 3.7 GeV^ are shown in Figure |13|, together with the 
NRQCD predictions. The colour-singlet channel underestimates the data by a factor ~ 
2-3, but there still is a considerable amount of uncertainty in the normalization from the 
value of the charm-quark mass, the wave-function at the origin, higher-order corrections 
as well as the choice of parton distribution functions and renormalization/factorization 
scale. The normalization of the theoretical prediction is significantly improved when 
adding colour-octet contributions; the differential distributions, however, seem to favour 
the shape predicted by the colour-singlet cross section. More theoretical work and 
more experimental information extending to larger Q^ values are needed before firm 
conclusions can be drawn. 

Various other observables and final states can be accessed with higher-statistics 
data, like photoproduction of tfj' and T particles, associated J /ip+'-^ production [ pB| , ^p^ 
as well as fragmentation contributions at large px [l68| , |59| . The future analyses at HERA 
will offer unique possibilities to test the NRQCD approach and to study the mechanism 
of quarkonium production in general. 
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Figure 13. Differential cross sections for ep — > eJ/tpX with Ef^^ > 5 GeV, 
Q^ > 3.7 GeV^ and p|i > 3.4 GeV^ |Q compared to theoretical predictions 
within NRQCD ||. 
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